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Introduction 

- Phase Retrieval / NASA Projects 

- JWSTTRL-6 
GBT Data / Notes: 

- Data Format and Sampling 

- Ray Trace Model & Wavefront 

- Symmetry of GBT Data 

- Pupil and Fourier Model 

- Pupil Amplitude 
PR Simulations 

- Wavefront derived from GBT Data symmetry 

- Wavefront Sensing accuracy and Coherent 
Assumptions 

GBT Results 
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>logy Development 
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Paiomar 200” Tefescope Adaptive Optics System (PALAO) Calibration 

10/2005 

JWST Testbed Telescope (TBT; Ball Aerospace, Boulder, CO) 
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1 estbed 1 elescope (TB i) 

Flight traceable, 1/6 scale, 18 segment design 

Algorithm Performance requirements dictated by NASA’s TRL-6 

- Testbed provides functionally accurate simulation platform for developing 
deliverable WFSC algorithms and software, 

- Used to perform TRL-6 end to end testing, 

- a solution is a fine-phasing algorithm that incorporates feedback, 

- an adaptive diversity function, eliminates Branch Points, and Wrapping 



Phase Retrieval Approaches 



• Two main approaches commonly used: 


- Iterative Transform (IT A) 

• Gerchberg-Saxton 

• Misell-Gerchberg-Saxton 

• HDA (extends dynamic range) 


- Parametric (non-linear least squares model fitting) 

• Solve for aberration coefficients 

• Solve for point-point phase in the pupil 



For JWST - adopted a 
hybrid approach that 
incorporates features of 
both types of algorithms. 


min [Objective Function]; 

— js Ij ?D.v J a -)|f -0 


f e g., W. H. Southwell, “Wavefront Analyzer using a Maximum Likelihood Algorithm," 
’ J. Opt. Soc. Am. A3, 396*399 (1 977). 
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Earlier Work using ITA with Radio Antennas 
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1988 , D. Morris, et al., “Radio holography measurement of the 30-m 
millimeter radio telescope Astron. Astrophys., vol. 203, p. 399. 

1991 , D. MORRIS, et. al, ‘Experimental assessment of phase retrieval 
holography of a radiotelescope’, IEE Proc. H, 138, pp. 243-247 

1994 , A. Greve, D. Morris, et. al., "Astigmatism in Reflector Antennas: 
Measurement and Correction," IEEE Trans ANTENNAS & Prop VOL. 

42, NO. 9 

1996 , D. Morris, Simulated Annealing Applied to the Misell algorithm for 
phase retrieval, IEE Proc - Microw Antennas Prop , Vol 143, No 4, 
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Notes / Understanding of GET Data 


• Consists of two feeds (pixels), two polarizations, 

• Separated by 58 arc-seconds, 

• Output of receivers is differenced to minimize the effect of sky- 
brightness variations. 

• Effective response of the telescope is modeled as the real beam 
convolved by two delta functions separated by 58” in the 
azimuth direction 

- aberrations due to both of the feeds being off (and on opposite sides of) 
the optical axis are negligible? 

- if this is not negligible, then a “single beam convolved by two delta 
functions” assumption may not be valid. 
















v = 43.1 GHz; A, = 6.96 mm 


Azimuth direction (x), approximately 350 samples/per scan line. 

Sampling in Azimuth = 3600*(180/pi)*(dx(251)-dx(250)) = 2.42 Arcsec / pixel = px = 1.1732e-05 radians 
Qx = X/(D*px) = 6.96e-3 / (100*1. 1732e-05) = 5.9325 

Elevation direction (y), 1 7 scan lines 

Sampling in Elevation = 3600*(180/pi)*(dy(5600)-dy(250))/17 = 7.5 Arcsec / pixel = py = 3.6361e-05 radians 
Qy = X /(D*px) = 6.96e-3 / (100*3.6361e-05) = 1 .9141 

Nyquist sampling is 7.2 arcsec / sample, Q = 2; Under-sampled by 0.96 in Elevation : over-sampled by 2.97 in 
Azimuth 







NRAO GBT Aperture 



• Panels are arranged in such a way that rings are concentric 
with a parent parabola. 


• Zemax design: GBT is 
setup as a single off-axis 
section of the parent 
parabola. 
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GB I F ouner Mode 



Note that Translation Shift of Fourier Transform produces a phase factor: 

3 [g(x-x 0 )] = G((O)e-‘ 2 ’‘“ x ‘ 



relative amplitude 


Pupil Illumination - 1 

Cool Link: edge taper in radio astronomy (Che 

http ://www .alma.nrao.edu/memos/htmi-memos 

A(r ) = 1, (uniform) 

= exp[-a (r / r 0 ) 2 ], ( tapered Gaussian ) 
a = (T e l 20)lnl0, (edge taper factor) 



meters 


3 / Mangum): 

frnal 97/memol 97.htm! 


T e = edge taper in dB 


----- A, 

amplitude variation: 

AA dB = 101og[(A 1 /A 2 ) 2 ] = 201og(l/0.18) 
= 15 
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relative amplitude 


Pupil 111 umination - II 


A(r) = 1, ( uniform ) 

yy r\ ry 

= exp[-r /2 <7 ], (tapered Gaussian) 
with a - 0.3, from PTCSPN47.pdf 

...the aperture plane amplitude distribution, that is, the illumination of the primary surface. This was 
approximated as a well-centered and circular Gaussian with a width (in radius-normalized units) defined 
by o = 0.3 . which corresponds to 1 4.5 dB of illumination taper at the edge of the dish . . . 

(7 = 0.55 

amplitude variation: 

A A dB = 1 0 log[(A , / A 2 f ] = 20 log(l / 0. 1 8 ) 
= 15 dB 



meters 
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Challenge for ITA Phase Retrieval 

Two incoherently subtracted irradiance values appear in the GBT data. 
Data collection process, I = I { ~ 1 2 

For the ITA approach to work, these irradiance values should be the 
result of one FFT. 

So make the approximation that: 


Coherent Approximation for Incoherent Data: 


3{A,(-0, HA. (+<?,)} =3 {A L (-e,)} + 3{A„(+0,)} 


or simply I ~ I L + L R 
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mage on 1 -side of focus: Dual aperture model 


Pupil Amplitude 


Results 





Co m ment on GBT Beam S 


Model 


Data 



Let PR solve for 
the wavefront 
that is consistent 
with the data: 



Incoherent Mode 











dual wavefront 


Incoherent Data 


Recovered. 


Pupil Amplitude 


Results 
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Lett Aperture 














S ummary 


• In principle, coherent IT A PR may work on 
incoherent GBT data, 

• Errors increases as beam tilt decreases 



